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Abstract. We present a new analysis of the minimum mass for star formation, based on opacity-limited fragmenta- 
tion. Our analysis differs from the standard one, which considers hierarchical fragmentation of a 3-D medium, and yields 
^mm ~ 1.007 to 0.010 M for Population I star formation. Instead we analyse the more realistic situation in which there is 
one-shot fragmentation of a shock-compressed layer, of the sort which arises in turbulent star-forming clouds. In this situation, 
M MIN can be smaller than 0.003M G . Our analysis is more stringent than the standard one in that (a) it requires fragments to have 
condensation timescales shorter than all competing mass scales, and (b) it takes into acount that a fragment grows by accretion 
whilst it is condensing out, and has to radiate away the energy dissipated in the associated accretion shock. It also accords with 
the recent detection, in young star clusters, of free-floating star-like objects having masses as low as 0.003 M Q . 
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1. Introduction 

3-dimensional opacity-limited hierarchical fragmentation 
yields a minimum stellar mass M mm ~ 0.007 to 0.010 M Q . 
This result was first obtained to order of magnitude by Rees 
(1976), but the range quoted above is the result of more 
detailed treatments of the thermodynamics of fragmenta- 
tion by Low & Lynden-Bell (1976), Silk (1977) and Boss 
(1988). However, there are problems with the process of 
3-dimensional hierarchical fragmentation, and this estimate of 
M MIN3 almost certainly needs to be revised upwards (see 
Moreover, it is now generally accepted (e.g. Elmegreen 2000, 
Pringle et al. 2001, Hartmann et al. 2001, Elmegreen 2002, 
Padoan & Nordlund 2002, Vazquez-Semadeni et al. 2003) 
that star-forming clouds are transient turbulent entities, and 
that star formation occurs almost as soon as a cloud forms, 
within a few crossing times. In this picture, prestellar cores 
are created where turbulent elements collide with sufficient 
ram-pressure for the resulting shock-compressed layer to be 
gravitationally unstable. The process of fragmentation is then 
quite different from that envisaged in 3-dimensional hierar- 
chical fragmentation. In particular, (i) the fragmentation of a 
shock-compressed layer can, and usually does, proceed whilst 
the layer is still accumulating (and therefore it is confined 
by ram pressure); (ii) the convergent motions leading to the 
formation of a fragment are initially concentrated in the plane 
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of the shocked layer, i.e. fragmentation of a shock-compressed 
layer is essentially two-dimensional. 

In this paper, we present a model to describe the 
2-dimensional opacity-limited fragmentation of a shock- 
compressed layer and use this model to determine the min- 
imum mass of the fragments that can condense out of such 
a layer, M MN2 . Section previews the standard paradigm of 3- 
dimensional hierarchical fragmentation, in which a cloud frag- 
ments into smaller and smaller subclouds until radiative cool- 
ing can no longer match the compressional heating rate and the 
the subclouds become approximately adiabatic. In Section[3]we 
highlight the differences between hierarchical 3-dimensional 
fragmentation and the fragmentation of a 2-dimensional layer. 
We derive equations to describe how a shock-compressed layer 
is formed by the collision of two streams of gas, and how the 
layer fragments whilst it is still forming. In Section 0] we use 
this model to follow the evolution of a fragment and hence to 
estimate the minimum mass for a fragment. In Section [5] we 
briefly summarize our results. 

Observational and numerical evidence for the formation 
and fragmentation of shock-compressed layers, in regions 
of imminent or ongoing star formation, is hard to identify. 
This contrasts with the substantial evidence for the formation 
and fragmentation of filaments in such regions. Nonetheless, 
we believe that the formation and fragmentation of shock- 
compressed layers must be the fundamental process in turbu- 
lent star formation regions. First (from a theoretical perspec- 
tive), to create a filament ab initio would require a very con- 
trived flow geometry. Second (again from a theoretical per- 
spective), the fragmentation of a layer normally proceeds via 
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the formation of a network of filaments (e.g. Turner et al. 1995, 
Whitworth et al. 1995, Bhattal et al. 1998); this is because of 
the well established theorem (e.g. Zeldovich 1979) which as- 
serts that if gravitational contraction gets ahead in one dimen- 
sion (due either to a greater initial perturbation, or less resis- 
tance), it tends to get further ahead. Third (from an observa- 
tional perspective), layers are much harder to identify than fil- 
aments, unless they are seen edge-on, in which case they look 
like filaments. 

Moreover, free-floating objects with masses as low as 
0.003 M Q have now been detected in young star clusters, like 
cr Orionis (Zapatero Osorio et al. 2002), and it beholds us to 
ask how they might have formed. They cannot form by three- 
dimensional fragmentation, but it appears that they could form 
by layer fragmentation, as analyzed here. Other possible for- 
mation mechanisms include the photo-erosion of pre-existing 
massive cores (Hester 1997, Whitworth & Zinnecker 2004) and 
impulsive interactions between discs in dense clusters (Boffin 
et al. 1998, Watkins et al. 1998a,b). 

2. Three-dimensional hierarchical fragmentation 

Theories of star formation have traditionally been based on 
the idea of opacity-limited hierarchical fragmentation (Hoyle 
1953). A massive, 3-dimensional, Jeans-unstable cloud starts 
to contract and, as long as the isothermal sound speed, a, re- 
mains constant, the Jeans mass 
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decreases with increasing density, p. Thus, once the original 
cloud has contracted sufficiently, the Jeans mass M JEANS3 is re- 
duced and the cloud can fragment into subclouds. Once these 
subclouds have contracted sufficiently, they themselves can 
fragment into still smaller 'sub-subclouds', and so on. The pro- 
cess can continue recursively, breaking the original cloud up 
into ever smaller fragments, as long as the isothermal sound 
speed remains approximately constant, i.e. as long as the con- 
tracting gas is able to radiate away, immediately, the PdV work 
being done on it by compression. Once the smallest fragments 
become so dense and opaque to their own cooling radiation that 
they can no longer radiate away the PdV work immediately, the 
sound speed starts to increase, the Jeans mass decreases no fur- 
ther, and fragmentation therefore stops. 

The minimum mass for star formation can be found from a 
general analysis of this process, as shown by Rees (1976); for 
simplicity we neglect all purely numerical factors. A fragment 
of radius R, which has just become Jeans unstable and started 
to condense out, contracts at speed dR/dt ~ -a. Hence the 
heating rate due to PdV work is 
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where we have substituted P = pa 2 and dV/dt ~ R 2 dR/dt. If 
the fragment is to remain approximately isothermal, this must 
be less than the maximum possible radiative cooling rate, 



where cr SE is the Stefan Boltzman constant and T is the temper- 
ature. This maximum cooling rate will only be realized if the 
fragment cools like a blackbody, i.e. if it is optically thick but 
only just; detailed calculations (e.g. Low & Lynden-Bell, 1976) 
indicate that this is usually the case. To ensure that C > 'H, i.e. 
that cooling is sufficently rapid to keep a constant, we require 
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or equivalently 
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Here we have substituted a = {k B T /m) 1 ^ 2 , where k B is 
Boltzmann's constant and m is the mean gas-particle mass; 
°"sb ~ k^lc 2 ^, where c is the speed of light and h is Planck's 
constant; and ot PL anck = [ch/G] 1 ^ 2 is the Planck mass 

However, there are problems with the notion of hierarchical 
fragmentation. Observationally, there is no clear evidence for it 
occuring, and from a theoretical viewpoint, it does not seem to 
work because the time-scale on which a subcloud condenses 
out is always longer than the time-scale on which its parent 
cloud is contracting. 

Specifically, the condensation timescale in 3-D is given by 

-,2m - 1/2 



1 - 



M 
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where f FF3 = {3n/32Gp) 1 ^ 2 is the freefall time in 3D. Thus a 
parent cloud is always closer to freefall collapse than its sub- 
clouds. If we assume that the whole of the parent cloud breaks 
up into subclouds, so that at their inception the subclouds are 
touching, then a subcloud has insufficient time to establish it- 
self as a distinct entity before it gets merged with neighbouring 
subclouds by the overall contraction of the parent cloud. For 
example, a subcloud with M ^ 4M JEANS3 condenses out on a 
timescale ~ 1.3 f FF3 (and smaller subclouds condense out even 
more slowly), whereas the parent cloud contracts on a time- 
scale ~ f FF3 . 

If, instead, we assume that at their inception the subclouds 
are not touching, we must consider whether a subcloud con- 
tinues to grow by accreting material from its surroundings as 
it condenses out. If we take at face value the Bondi formula 
for the rate of spherically symmetric accretion onto a mass M, 
from a background medium having density p and isothermal 
sound speed a, 



dM 
dt 



e V2 jrG 2 pM 2 



(8) 



C 



R 2 cr„ 



(3) 



1 Equation ^6) shows that M MIN3 is much smaller than the 
Chandrasekhar mass by virtue of (i) the mean gas-particle mass m 
being much larger than the mean mass per electron (2ra pROTON /(l +X), 
where X is the fractional abundance of hydrogen by mass); and (ii) the 
sound speed, a, being much less than the speed of light, c. 
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(Bondi 1952), we can rewrite this equation in the form 



din(M) 
d[t/t m ] 
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(9) 



Thus, even though the Bondi scenario is not strictly applicable, 
the implication is that the subcloud mass will increase by a very 
large factor whilst it is condensing out (e.g. a fragment with 
initial mass 4M JEANS , notionally increases its mass by a factor of 
42 x 4 x 1 .3 ^ 220). We must conclude that the initial fragment 
mass will be a significant underestimate of the final fragment 
mass. 

3. Fragmentation of a two-dimensional layer 

There is an alternative paradigm to hierarchical fragmentation, 
which overcomes these problems and agrees better with obser- 
vation. In this paradigm (e.g. Larson 1981, Elmegreen 2000, 
Pringle, Allen & Lubow 2001, Hartmann, Ballesteros-Paredes, 
& Bergin 2001, Padoan & Nordlund 2002, Elmegreen 2002, 
Vazquez-Semadeni, Ballesteros-Paredes & Klessen 2003, Mac 
Low & Klessen 2004), giant molecular clouds are relatively 
short-lived objects which form and dissolve on a dynamical 
timescale and have a turbulent and inhomogeneous internal 
structure. The substructure within a molecular cloud is there- 
fore highly transient, with clumps forming, dispersing, and re- 
forming on a dynamical timescale, without necessarily spawn- 
ing new stars. Only occasionally will a particularly dense, mas- 
sive and strongly converging shock lead to the formation of a 
prestellar core. In this scenario, both core formation and core 
collapse are triggered by the collision of two turbulent ele- 
ments of gas and the formation of a gravitationally unstable, 
shock-compressed layer. This is a one-step, two-dimensional 
fragmentation process: 'one step' because the process does not 
repeat itself hierarchically; 'two-dimensional' because, when 
the layer fragments, the wavelengths of the most rapidly grow- 
ing fragments are much larger than the thickness of the layer 
(Whitworth et al. 1994a,b). 

3.1 . Linear fragmentation of a static layer in 
plane-parallel symmetry 

2-dimensional fragmentation of a static layer (e.g. Larson 
1985) is fundamentally different from 3-dimensional fragmen- 
tation. The Jeans mass is given by 
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where 2 is the surface density of the layer; and the timescale 
on which Jeans unstable fragments condense out of the layer is 
given by 
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Since f COND , has a minimum for M 



4M jeans2' fragments on 
this mass scale condense out faster than smaller fragments 
and faster than larger fragments. Therefore 4M JEANS2 is a pre- 
ferred mass scale for fragmentation. Fragmentation on this 



scale is likely to be permanent, because the fragments will 
not be merged by the overall contraction of the parent layer 
on larger scales. Given these differences between hierarchical 
fragmentation of a 3-dimensional cloud and one-step fragmen- 
tation of a 2-dimensional layer, we are interested in exploring 
how, and at what mass scale, opacity limits the fragmentation 
of a 2-dimensional layer. Is the minimum mass formed at the 
opacity limit significantly different from that estimated for 3- 
dimensional fragmentation? 

3.2. A shock-compressed layer in plane-parallel 
symmetry 

We consider two identical streams of gas with uniform density 
p and uniform isothermal sound speed a, which collide head- 
on at relative speed 2v to form a layer. We assume that a plane 
contact-discontinuity forms where the two streams meet and 
we fix our main coordinate frame in this contact discontinuity. 
The two streams approach the contact discontinuity with veloc- 
ities ±v, and a plane-parallel layer forms symmetrically about 
the discontinuity, bounded by two accretion shocks, as illus- 
trated in Figure 13.21 For simplicity we assume that radiative 
cooling in the shocked gas is so efficient that the gas rapidly 
cools back to its pre-shock temperature. If the full thickness of 
the layer is 2Z(f), then the gas flows into the accretion shock at 
speed K m = v + Z, where Z = dZj dt, and out at speed m out = Z; 
note that u m and m out are measured relative to the shock-front. 
Applying the isothermal shock condition, 

M„, fl 

— = , (12) 



we obtain 

v + Z a 

a 2' 

whence 

[v 2 +4a 2 ] 1/2 -v 
Z = 



(13) 



(14) 



The post-shock density (i.e. the density in the shock com- 
pressed layer) is 



Pa 



= P" 



= P 



[v 2 +4a 2 ] 1/2 + v 
[v 2 +4a 2 ] 1 / 2 -v' 



(15) 



and so the surface density 2 of the layer grows linearly with 
time according to 



2(f) = 2p SHOCKED Zf = p{[v 2 +4a 2 ] 1/2 + v}f. 
As long as v » a, we can approximate 



Z =s — <K a , 
v 

2 

p V 



2(f) =s 2pvt, 



(16) 

(17) 

(18) 
(19) 



and so the sound-crossing time of the layer, Z/a, is much 
smaller than its growth time, Z/Z. This means that the layer has 
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sufficient time to relax and remain close to hydrostatic equilib- 
rium. We also note that, provided 



t < 



1 



2(Gp) 



1/2 



(20) 



GZ 2 «: p m0CKED a 2 , and thus self-gravity plays a negligible role 
compared with ram pressure in confining the layer in the direc- 
tion perpendicular to the contact discontinuity. Consequently, 
the plane-parallel hydrostatic equilibrium of the layer has a 
very flat density profile, with little density contrast between the 
centre and the edge. 



3.3. Linear fragmentation of the layer 

Now we consider a small circular patch of radius r on the layer, 
and consider whether it is able to condense out. The radial mo- 
tion of the patch is determined by a competition between the 
hydrostatic acceleration ~ a 2 jr (which promotes expansion), 
and the self-gravitational acceleration ~ nGI. (which promotes 
contraction). Therefore we can write 



r nGI, 

r 

and contraction (r < 0) requires 

a 2 

r > R 



nGI 



(21) 
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The time-scale for condensation of a Jeans-unstable fragment 
is given by 



1/2 



nGI 



,-1/2 



(23) 



so the fastest condensing fragment has initial radius and con- 
densation time-scale 



r =s 2R 

FASTEST JEANS2 



2a 2 
nGI' 



2a 
nGI 



(24) 



(25) 



As the layer piles up, I increases, and so the size of the 
fastest condensing fragment decreases, and its condensation 
time-scale decreases. Condensation into the non-linear regime 
starts only when ? FASXESX S= t, and substituting for 2 in equation 
J25b . from equation d!6i . this condition gives the fragmentation 
time, 
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1/2 
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The mean radius and mass of the fragments which condense 
out are therefore 



r (^FRAG ) 
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We note that 
2r 



2a 



2Z( tFRAG ) [( v 2 +4fl 2 ) 1/2 -v] 



(29) 



where the final expression obtains in the limit v » a. 
Therefore, at its inception, the mean diameter of a fragment is 
much larger than the thickness of the shock-compressed layer. 

3.4. Non-Linear fragmentation of the layer 

The linear analysis above does not take into account the ma- 
terial which continues to flow into the shock layer, after it has 
started to fragment. It also does not describe how the fragment 
develops as it condenses out. 

In order to follow the condensation of fragments into 
the non-linear regime, we model them as uniform density 
spheroids with radius r and half-height z. The excursions of 
an isolated oblate spheroid of mass m, subjected to external 
pressure P EXT , are governed by the equations 
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In equations J30l > and ( 13 U . the first term on the righthand 
side represents self-gravity, the second term represents ex- 
ternal pressure, and the third term represents internal pres- 
sure. Similar expressions describe the excursions of a prolate 
spheroid. 

For an oblate spheroidal fragment condensing out of a 
shock-compressed layer, while the layer continues to accrete 
matter from the converging flows which are forming it, we must 
modify these equations to take into account two facts, (i) The 
external pressure is the ram pressure of the inflowing gas. (ii) 
The spheroid grows in mass as a consequence of the inflowing 
gas. Knowing z, which is the speed with which the gas already 
in the spheroid expands, we can calculate how fast, y, the shock 
front has to advance towards the inflowing gas, in order to de- 
celerate the new inflowing gas. For the velocities relative to the 
shock, we have « IN = v + y and u om = y - z, so the isothermal 
shock condition gives 



v + y a 



a y-z 



[(v-z) 2 + 4a 2 ] 1 ' 2 - (v - z) 



y 



(32) 



(33) 



The external pressure acting on the spheroid in the z-direction 
is therefore 



^EXT. ; = P° 



2 v + y 

y-z 
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Fig. 1. Gas streams having density p and sound speed a collide to create a shock-compressed layer. As gas continues to flow into 
the layer from each side at speed ±v, a spheroidal fragment of initial radius and initial height z !NIT = Z begins to condense 
out. 



[(v-z) 2 + 4fl 2 ]' /2 + (v + z) 
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The external pressure acting on the spheroid in the r-direction 
is the same as in the unperturbed layer, 

|l/2 



pa 



v 2 + 4a 2 ] 



+ V 



(35) 



[v 2 +4a 2 ] 1/2 -v' 

The rate of increase of the mass in the fragment, due to 
accretion from the continuing inflow, is 



2 7r r p (v + y) 



xr 2 p {[(v-z) 2 +4a 2 ] I/2 + (v + z) 



(36) 



and we must add a term to equation (I30> to represent the inertial 
drag of this accreted material, i.e. 

3Gm (rcos^iz/r) (z/r) 



(r 2 - z 2 ? 12 (r 2 - z 2 ) 

20nP rYT rz 5 a 2 mr 

r^ 2 — + . (37) 

5m r m 

The PdV heating rate for the material already in the frag- 
ment is 



<H = -P— = -ma 2 
dt 



2r z 
h - 

r z 



(38) 



The heating rate due to dissipation of the z-kinetic energy of 
the matter accreting onto the fragment from the inflow, less the 
r-kinetic energy used to accelerate the accreted matter laterally, 
is 



m (v + z) 



m r 



(39) 



The maximum (i.e. blackbody) cooling rate from the two 
sides of the fragment is 

2 2 tt 6 m 4 a & r 2 



C(T) = 2;rrV SB T = 



(40) 



15 c 2 h 3 ' 

where h is Planck's constant and ra is the mean gas-particle 
mass, as before. 



3.5. Survival criterion 

We can now follow the evolution of a collapsing fragment. The 
opacity limit for the proto-fragment is reached as soon as "H + 
T> <: C(T). In reality, the radiative cooling is more complicated 
than equation (I40> . The gas flowing into the accretion shocks 
which define the top and bottom of the fragment (see Fig. 1) is 
initially heated to high temperature (T sli0CKED ~ 3mv 2 /16k B ) in 
the shock front, and then rapidly cools behind the shock front. 
Provided we impose the condition pv 3 /2 <k o" sb ^ , or 



p « p c 



t3 5-4 5 

3 7r m v 
2 14 5c 2 /! 3 



KT 12 gcnr 3 



kms 



(41) 



then the immediate post-shock cooling (IPSC) radiation is ei- 
ther optically thin continuum radiation (e.g. from dust), or op- 
tically thick cooling in a few specific molecular lines which 
occupy a small total bandwidth. In either case, half the IPCS 
radiation is radiated away from the fragment, and half towards 
it. In the most extreme case, all of this latter half is absorbed by 
the fragment, and then has to be re-radiated at temperature T, 
along with the internal energy delivered by PdV compression. 
Thus, by requiring the fragment to radiate at temperature T all 
the energy dissipated in the accretion shock, we are making a 
rather conservative assumption. 

To ensure that a proto-fragment forms a distinct condensa- 
tion, i.e. that it continues to undergo contraction after the opac- 
ity limit has been reached, we set the following two survival 
conditions at the opacity limit: 

1 . the proto-fragment must be contracting in both the radial 
and vertical dimensions, i.e. 



< 0, 



LIMIT 

z,™,t < 0; 



(42) 
(43) 



2. the final extent of the proto-fragment, in both dimensions, 
must be less than or equal to half its initial radial size, i.e. 



< 5 r 

LIMIT yj.^ i INrT ■ 

< 5 r 

LIMIT yJ " J ' INIT * 



(44) 
(45) 
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Fig. 2. A log/log plot of the (p, v) plane. The dots mark combinations of pre-shock density, p, and collision speed, v, for which 
(assuming a post-shock effective sound speed a = 0.2kms~', corresponding to molecular gas at 10K) the fastest-growing 
fragment mass is less than 0.005M o . The irregularities in the boundary of this region have to do with the tendency of fragments 
with mass near the fastest growing one to undergo pulsations - see 34.21 and 34.31 The solid line is the locus (equation ( 14 11 ) 
below which p must fall if the immediate post-shock cooling radiation is to be negligible (as we assume). 



For given p, v and a (i.e. for a given colliding flow), we must 
pick f INIT and r 1N1T , where f INIT is the time at which the fragment 
under consideration starts to condense out, and r INIT is its initial 
radius. Thus the initial conditions for a trial fragment are 



t = f INIT , (46) 

'"(^INIt) — '"iNIT ' (47) 

Kt mn ) = 0, (48) 

z(f INIT ) = Z(f INIT ), (49) 

zit^) = Z(f INIT ). (50) 



We can then use equations J37L ( 13 11 . i34\ . J35i and d36l > to 
follow the development of the fragment. Then we can use 
equations J38l >. d39l and J40i to evaluate H, D, C, and check 
whether the fragment has reached the opacity limit (H + T> <: 
C). Finally, once the limit is reached, we can check the survival 
conditions, equations i42i . t43l . J44l > and ( 145 \ . 

In order to reduce the parameter space we fix a = 
0.2 kms -1 , correponding to molecular gas at 10 K. The equa- 
tions are integrated using a fourth-order Runge-Kutta scheme. 

In reality the gas in the two colliding streams will not have 
uniform density, and so p should be interpreted as the mean 
pre-shock density. The presence of inhomogeneities in the pre- 
shock gas will promote fragmentation of the resulting layer, by 



creating the seed perturbations from which condensations sub- 
sequently grow. It will also tend to yield condensations with 
finite angular momentum. However, the lowest-mass conden- 
sations will be those with low angular momentum, and so we 
ignore angular momentum in this exploratory analysis. 

4. Results and discussion 

4.1. General considerations 

A condensing fragment posseses two orthogonal dimensions r 
and z. Moreover, (a) r and z in general have different initial 
values (r(f INIT ) and z(f INIT )). (b) Their excursions are driven by 
different external pressures (see equations \3A\ and \35Y ). (c) 
The coupling between them (which occurs through both inter- 
nal pressure and self-gravity) is non-linear. Therefore a frag- 
ment can evolve in quite complicated ways. In particular, small 
fragments which are not initially very unstable tend to oscillate 
until they accrete sufficient mass from the continuing inflow to 
become unstable against monotonic condensation. The oscil- 
lations in r and z normally have unrelated periods, and so the 
model has to be able to treat both oblate and prolate spheroids. 
However, once a fragment starts to condense out it is normally 
oblate, due to the extra ram pressure exerted in the z dimension 
by the inflow. 
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c) 



e) 



g) 



1.5x10"° 2x10"- 
t(Myrs) 




b) 



d) 




h) 



l.5x10" J 2x10 
t(Myrs) 





i) 

Fig. 3. Evolution of the minimum-mass fragment: (a) radius, r(f) in AU; (b) height, z(t) in AU; (c) radial velocity, r(f) in units 
of a = 0.2kms _1 ; (d) vertical velocity of material already in the fragment, z(t) in units of a; (e) vertical velocity of shock 
front bounding the fragment, y(t) in units of a; (f) mass of fragment, m(t) in M JUPITER ; (g) accretion rate, dm/dt in M yr _1 ; (h) 
logarithmic fragment density, £og(p FRAG ) in g crrT 3 ; (i) logarithmic ratio of heating to cooling, £og(['H + D]/C). 
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4.2. Oscillations 

If we fix p and v, and vary f INIT and r INIT , we find that the lowest- 
mass condensation, m.,,.,(p, v), forms when f__.„ and 

1 MIN ' ' 1 INIT FRAG 

r INIT =a 0.5r FRAG . Fragments with r INIT ^ 0.5 r FRAG tend to un- 
dergo radial oscillations before accumulating sufficient mass to 
condense out. Fragments with r 1N1T > 0.5 r FRAG tend to condense 
out monotonically. The lowest-mass fragments are those which 
start with just enough mass to condense out monotonically, i.e. 
^init ~ 0.5t*p RAG . Since 0.5 f" FRAG 

- ^jeans2' this finding agrees 
with the predictions of the linear stability analysis in 33.31 

4.3. The optimum combination of p and v for 
producing low-mass condensations 

Suppose now that only p is fixed, and v is varied. For each v 
we can determine the lowest-mass condensation m MIN (p, v), by 
varying f INIT and r INIT as described above. What we find is that 
as v is increased at fixed p, m MIN (p, v) at first decreases, reaches 
a minimum, and then increases. Thus for any value of p there 
is an optimum value of v for spawning low-mass condensa- 
tions. Fig. 2 shows the (p, v) plane and marks discrete points 
for which m wm {p, v) is less than 0.005 M = 5M JUPITER . We see 
that there is a significant region of (p, v) space in which con- 
densations with mass below 0.005 M can form. 

This region has an irregular boundary because proto- 
fragments can undergo out-of-phase pulsations in the r- and 
z-dimensions. If a proto-fragment undergoes a pulsation in the 
r-dimension and expands for a while (r > 0), it experiences a 
prolonged period of growth due to accretion onto its relatively 
large cross-section (nr 2 ), and hence its mass increases rapidly. 
Since the onset, amplitude and phase of such pulsations is criti- 
cally dependent on initial conditions, small changes in p and/or 
v can result in quite large changes in m MIN . 

The full line on Fig. 2 shows the locus p = p CRIT = 
10~ 12 gcm~ 3 [v/kms 1 ] 5 . Hence the condition for the IPSC ra- 
diation to be neglected (equation (14 1» requires p to be well 
below this line. Evidently this condition is always easily ful- 
filled. 



quickly resumes its radial contraction. Soon after this the opac- 
ity limit is reached, and we must presume that the fragment 
then switches to Kelvin-Helmholtz contraction. At this stage 
the fragment is roughly spherical. It is still accreting, at a rate 
~ 2 x 10~ 4 M JUPITER yr _1 , and it has a dynamical timescale of 
~ lOOOyr, so we might expect its mass to increase by a further 
few tenths of a Jupiter mass during the subsequent condensa- 
tion. 

5. Conclusions 

In 3 -dimensional hierarchical fragmentation, the minimum 
mass is estimated to be M mm ~ 7 to 10M JUPITER (e.g. 7M JUPITER , 
Low & Lynden-Bell 1976, 10M JUPITER , Silk 1977; 10M IUPITER , 
Boss 1988). Moreover, these estimates probably need to be re- 
vised upwards, to take account of ongoing accretion and merg- 
ing during concensation (see |2j. 

However, surveys of young clusters are beginning to find 
objects with masses estimated to be as low as 3 M JUPITER (e.g. 
Zapatero Osorio et al. 2002). 

We find that, in the 2-dimensional fragmentation of a shock 
compressed layer, the minimum mass is significantly smaller 
than 10M JUPITER . For certain rather specific shock parameters, 
fragments with mass < 3 M JUPITER can condense out - specifi- 
cally M sml — 2.6 M JUPITER - and for a wide range of shock pa- 
rameters fragments with ^ 5M JUPITER can condense out. Since (i) 
layer fragmentation is probably a more realistic model for frag- 
mentation in highly turbulent, transient star-forming clouds, 
(ii) layer fragmentation is more likely to be permanent than 
3-dimensional fragmentation (because the fragments in a layer 
condense out faster than the layer as a whole), and (iii) our anal- 
ysis takes proper account of ongoing accretion during fragmen- 
tation, we conclude that free-floating planetary-mass objects 
can form by the fragmentation of a shock-compressed layer, 
in the same way as stars and brown dwarves. 

Acknowledgements. DFAB thanks PPARC for the support of a post- 
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4.4. The minimum mass 

The lowest-mass condensation of all is formed when p st 9.6 x 
10- I6 gctrr 3 , v 5.1 a 1.02 km s -1 , ? INIT = f FRAG lOOOyr, 
and r INIT = 0.62 r FRAG =s 25 AU; the initial half-height of the 
fragment is z(? INIT ) - 8 AU, and the initial mass of the fragment 
is m(t [Kn ) - 1.5M JUPITER . By the time the fragment has reached 
the opacity limit, its radius is ~ 12 AU, and it is roughly spheri- 
cal, but slightly oblate. 1600yr have passed since its inception, 
and its mass has increased by 1.1M JUPITER to 2.6M JUPITER . Thus 

^MIN2 — ^JUPITER" 

Fig. 3 shows the development of this fragment. For the first 
~ 1200 yr, the fragment contracts in the radial dimension, but 
expands in the vertical dimension, due to accretion of extra 
material from the inflow into the layer. Around 1200 years af- 
ter its inception the fragment starts to contract in the vertical 
dimension. At the same time there is a very small bounce in 
the radial dimension, but this is short-lived, and the fragment 
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